Quantitative characterization of non-DLVO factors in the aggregation of black soil colloids

The variable role and fate of soil colloids under different environmental conditions are derived from their dispersion and aggregation properties. In this work, dynamic and static light scattering were used to characterize the original size, aggregation kinetics of natural black soil colloids (BSCs) and structural features of aggregates in electrolytes with different cations (K+, Mg2+, Ca2+), respectively. For these three cations, the aggregation kinetics followed the trend of Ca2+ > Mg2+ > K+ and the critical coagulation concentration (CCC) followed the sequence: K+ (134.30 mmol L−1) > Mg2+ (13.27 mmol L−1) > Ca2+ (4.19 mmol L−1). The results indicated that the aggregation behavior in different valence cation systems followed the classical Derjaguin-Landau-Verwey-Overbeek (DLVO) model qualitatively. However, the quantitative differences of CCC suggest the existence of ion-specific effects. The effective ionic charge coefficient 1.31, 2.20, and 2.78 of K+, Mg2+ and Ca2+ were proposed to consider of all the non-DLVO factors, which were obtained by forming a relationship based on mathematic between the electrostatic repulsion and the van der Waals attractive interaction at the CCC. The non-classical polarization of cations in a strong soil electric field is a primary mechanism of cation effects on soil colloid interactions, causing the difference in colloid interaction energy and further affecting soil colloid aggregation. This result is crucial for enriching the theory of charged colloidal interactions.

www.nature.com/scientificreports/ theory 23,24 . Nevertheless, several studies have also shown that critical coagulation concentration (CCC) even in the same valence ions inducing the aggregation of nanoparticles can be able to show different 25,26 , which is referred to as the specific ion effect or the Hofmeister effect. The current theories state that for a charged colloidal particle system, strong cationic polarization in cation-surface interactions strongly influences the surface properties of the colloid and the interaction energy between particles and thus, govern the stability of the colloid suspension 11,27 . The strong polarization increases with the increase in electric field strength or decreases with the external cation concentration 27 . Thus, polarization of ions at the interface may also be considered as one of the non-DLVO factors.
In general, it is customary to select one or two separate mineral or organic soil components as model materials. However, the mature soil is poly functional and contains compounds with diverse composition and structure. Therefore, it is appropriate to conduct research on the aggregation behavior using real soil particles. The black soil in northeast China mostly has both permanently and variable charge on their surface because of a high content of mineral and organic components. The presence of natural organic matter, with a soft multibranched structure and high molecular weight and surface charge number, can both significantly inhibit the aggregation potential of mineral colloids/nanoparticles 28,29 , and also act as a binding substance promoting the formation of soil aggregates in natural soils 30,31 . The electric field strength at the surface of natural soil particles can reach 4 × 10 8-9 V m −1 . 32,33 . Some pertinent research questions which need to be addressed include the followings: What's the aggregation response of the black soil particles with higher organic content in strong soil electric field? How strong the non-DLVO factors affect soil colloid dispersion and soil aggregate structure?
The purposes of this study were to: (1) clarify the specific ion effects of K + , Mg 2+ or Ca 2+ on natural soil colloid aggregate stability; (2) quantitatively characterize the effect of non-DLVO factors on particle interaction while considering the specific ion effects and try to explore their main sources. The results can strengthen the understanding of the effects of ion species on the colloid aggregation process.

Materials and methods
Preparation and characterization of black soil colloids. The soil colloids were extracted from the black soil (Luvic Phaeozem, FAO soil classification) in Harbin (45° 19′ N, 126° 18′ E), Heilongjiang Province, Northeast of China. The black soil developed from a loess-like parent material and Maize (Zea mays L.) is the principal crop of this region. The soil sample from the top 0-20 cm layer was collected and sieved through a 0.25 mm mesh after air dried. Then the black soil colloids (BSCs) were prepared following the procedures 14,34 : first, 50.0 g black soil was transferred into a 500-mL beaker containing 500 cm 3 of deionized water. After thoroughly stirring the mixture, the pH of the suspension was adjusted to 7.5 ± 0.1 by 0.5 mol L −1 KOH solution. Second, the soil suspension was intensively sonicated (Scientz-IID, Ningbo, China) for 15 min and then diluted to 5000 mL with ultrapure water. Third, the < 200 nm (hydrodynamic diameter) fraction was extracted following the static sedimentation method 14 , and electrodialysis was used to remove free ions from the solution. A particle concentration of ~ 0.8185 g L −1 was determined by the oven-dry method. The specific area and CEC of BSCs determined by combined measurements are 180 m 2 g −1 and 24.2 cmol (+) kg −1 , respectively 35 . Which are now ready for dynamic light scattering (DLS) experiments. Moreover, < 2 μm soil colloid suspension was extracted also for measurement of the zeta potential.
Dynamic and static light scattering measurements. DLS was used to measure the hydrodynamic diameter of particles. A BI-200SM multiangle laser light scattering instrument (Brookhaven, USA) was used for the DLS measurements. All the measurements were carried out at wavelength of 532 nm and scattering angle of 90°. The scattering instrument was switched on to warm up for at least 25 min before an experiment. A total of 2 cm 3 of the BSC suspension was add into the scattering bottle. In order to attain the electrolyte concentration that experiment required and a total volume of 10 mL, moderate electrolyte solution and deionized water should be added to the bottle too. The BSC particle concentrations in the suspension were 0.1637 g L −1 . The average hydrodynamic diameter of the original BSCs measured by DLS was 192.0 nm (Fig. 1). Time-resolved DLS was used to measure the change in the number-weighted average hydrodynamic diameter of BSCs and their aggregates as a function of time in the presence of different cations. The size of the BSC aggregate was recorded every 30 s for 60 min. All measurements were carried out at 25 ℃.
Determinations of the critical coagulation concentration (CCC) and the activation energy of aggregation. The stability of the colloidal suspension in the electrolyte solution is usually characterized by CCC 36 . Researchers established a method for determining the CCC of a polydisperse nonspherical colloid suspension by DLS measurements 36 . Then, the activation energy of particle aggregation was further measured 27 . In this study, we measured the total average aggregation rates (TAA) at different electrolyte concentrations, and then estimated the CCC. The TAA rate was defined as per Eq. (1): where, ṽ T (f 0 ) (nm min −1 ) is the rate of change of TAA in the time range from t = 0 to t = t 0 , f 0 (mmol L −1 ) is the electrolyte concentration, D(t) (nm) and D 0 (nm) represent the aggregate hydrodynamic diameter at t and t 0 , where, ΔE(f 0 ) (J mol −1 ) is the activation energy, T (K) is the absolute temperature, k (J K −1 ) is the Boltzmann constant, and ṽ T (CCC ) is the TAA rate at f 0 = CCC.
Electron microscopy and static light scattering characterize. The surface morphology and elemental composition of aggregates were carried out using field-emission scanning electron microscopy (FESEM) (Zeiss Ultra Plus). Specimens were sprayed with gold in a vacuum coating machine for 120 s, and examined at 2-30 kV. The analysis showed that the semi-quantitative elemental contents of BSCs were 24.74% of Si and 39.33% of O, followed by 13.26% of Al, 8.98% of Fe, 7.17% of C, 1.06% of Ca, 0.73% of N, and 0.23% of Mn. Predominant clay minerals were illite and montmorillonite. The static light scattering (SLS) measurement was used to characterize the structure of aggregates formed under different electrolyte conditions. The scattering index was obtained by measuring the change of the scattering light intensity with the scattering vector in the BSC aggregation process. When the scattering index is a constant, it is the fractal dimension (d f ) of the aggregate. The d f was used to characterize the structural compactness of the formed aggregates.
Zeta potential measurements. The zeta potential of soil particles was measured by a Zeta Plus (Brookhaven, USA). Negatively charged surfaces were observed for the black soil colloid with a zeta potential value of − 26.34 mV at pH = 7. The suspensions were prepared by adding 1 cm 3 of the ˂ 2 μm BSC suspension (1.255 g L −1 ) into 9 cm 3 of different species and concentration electrolyte solutions, namely, 30-350 mmol L −1 KCl and 0-10 mmol L −1 CaCl 2 and MgCl 2 under neutral pH conditions. Each measurement was performed 3-10 times at room temperature (25 ℃).
To estimate long-range van der Waals force, Eq. (9) is the value used by the aid 39 .
where, A (J) is the effective Hamaker constant. In this study, we use 6.99 × 10 -20 J as the Hamaker constant for black soil colloid 25 . From this equation, the P vdW (λ) depends only on the inter-particle distance and the Hamaker constant.
In simulating the DLVO force of BSC aggregation in K + system, the surface potential φ 0 of colloid particles can be obtained by Eqs. (6) and (7). The potential value is then brought into Eq. (5) to compute the middle point potential φ(λ/2). Then, this φ(λ/2) value is used to compute the P EDL (λ) at different electrolyte concentrations through Eq. (4). Finally, the net DLVO force is obtained by summing the P EDL (λ) and P vdW (λ) obtained by Eq. (9). While in Mg 2+ and Ca 2+ electrolyte systems, the DLVO force calculation process is the same as K + system, except that Eq. (8) is used to calculate φ 0 instead of Eqs. (6) and (7).
Thus, the energy barrier ΔW (kT) for particle aggregation is expressed in Eq. (10): 16: where, In the P(λ) ~ λ curve, the lower limit when P(λ) = 0 is defined as x (m), and the upper limit is defined as y (m). In addition, there is a potential barrier in the domain of λ = x to y. s (m 2 ) represents the average area passing through the repelling space in this interval from λ = x to y.

Results and discussion
Aggregation kinetics of BSCs in monovalent and divalent electrolytes. The average hydrodynamic diameter of BSC aggregates in the presence of KCl, CaCl 2 and MgCl 2 were measured at a range of electrolyte concentrations to explore the electrolyte-dependent aggregation heterogeneity (Fig. 2). According to the effective collision probability, two aggregation mechanisms for colloidal particle aggregation can be summarized. One mechanism, corresponding to fast aggregation, is considered as diffusion limited cluster aggregation (DLCA). The other, corresponding to slow aggregation, is the reaction limited cluster aggregation (RLCA). In the RLCA mechanism, the aggregate diameter increases with time are linear dependence, while it increases exponentially with the increase in time in the DLCA mechanism 25,40-42 . The aggregation of BSCs increased with the increase in electrolyte concentration because the cations neutralized the negative surface charge of BSCs, thereby reducing the repulsive force between particles and enhancing aggregation. Divalent cations (Ca 2+ and Mg 2+ ) could induce a rapid and an intense aggregation of BSCs than monovalent cations (K + ) even at far lower ionic concentrations. For example, the BSCs remained stable in K + systems at an ionic concentration of 10 mmol L −1 , while rapid aggregation occurred in Ca 2+ and Mg 2+ systems at the same ionic concentration; additionally, Ca 2+ had a stronger effect in comparison to that of Mg 2+ . Compared with monovalent cations, a much lower concentration of divalent cations induced a drastic BSC aggregation, which may be because high-valent cations can make the surface charge of BSCs less negative through the increased charge screening effect, which benefits aggregation 43 . Previous studies have found that the stronger the surface charge, the more significant the shielding effect of high-valent cations 18 . Contrary to the Schulze-Hardy rule, cations with the same valence can have a significantly different coagulation ability 1,44 . The Schulze-Hardy rule, emphasizing the theoretical coagulation value as determined by the (5)  26,27,34 . The data presented herein show that the BSC aggregation process was the strongest for Ca 2+ , followed by Mg 2+ , the lowest for K + . In summary, the dispersion stability of BSCs for the same electrolyte concentration were higher in KCl than that in MgCl 2 and CaCl 2 solutions. The effects of these cations on BSC aggregation are in the order Ca 2+ > Mg 2+ > K + , suggesting that there are strong specific ion effects in the BSC aggregation process.
FESEM images and d f of BSC aggregates with K + , Mg 2+ and Ca 2+ electrolytes. The different aggregation process not only affect the dispersion stability of soil colloids but also alter the structural characteristics and surface topography of aggregates. Thus, the d f was measured to characterize the BSC aggregate structure.
A relatively high d f indicates that the aggregate structure is dense and the porosity is small; a low d f indicates that the structure is loose, open and has a large porosity. The fractal dimension of aggregates under different concentrations of K + , Mg 2+ and Ca 2+ electrolytes were list in Fig. 3. Under extremely low electrolyte concentrations, the aggregation process is very weak. At this time, the aggregates size is small with insignificant fractal characteristics, and the fractal dimension is low. With the increase of electrolyte concentration, the aggregation process intensifies, and forming large aggregates with significant fractal characteristics. However, when the electrolyte concentration further increases, the repulsive barrier between particles disappears completely. At this time, the aggregation occurs immediately once the particles collide. Which is not conducive to form a regularly arranged fractal structure, so the fractal dimension of aggregates decreases. This result is consistent with the conclusion that RLCA contributes to the formation of aggregates with less porosity, compact structure and more prominent self-similar characteristics than DLCA 40,41 .
The data in Fig. 3 shows that the d f of BSC aggregates in Mg 2+ were all higher than that in Ca 2+ at the same ionic concentration. This trend suggests that ion-induced week coagulation is more likely to cause a compact aggregate structure. For black soil colloid aggregates, the reason for the dense soil structure is the dense aggregates, and a denser structure of Mg 2+ -aggregates means stronger repulsive forces among soil particles in soil solutions. These repulsions create opportunities for the particles to be bonded through multiple collisions, thus, finding more stable sites. Therefore, Ca 2+ -aggregates with lower fractal dimensions may correspond to a less stable soil structure. On the other hand, the fractal dimension of Ca 2+ -aggregates is low, which represents larger aggregate volume and stronger anti-sticking properties, viscosity resistance as the aggregates move in the aqueous medium of soil pores; therefore, aggregates with lower fractal dimension may travel slower through a soil pore 14 . www.nature.com/scientificreports/ The texture and shape of BSC aggregates, formed in different electrolyte solutions observed in FESEM and shown in Fig. 4, indicating that the aggregates are floc-like in Ca 2+ solution, schistose-like in Mg 2+ solution and block-like structure in K + solutions, with small cracks and variable degrees of roughness. Furthermore, the BSC aggregates induced by K + has a relatively smooth surface and clean edges. Comparatively, Mg 2+ and Ca 2+ -induced aggregates take the shape of sheet or floc and have a rough edge, with small bumps and variable degrees of roughness. It is evident that the structural changes and enhanced aggregation of BSCs are primarily due to the cation-surface interactions in a strong electric field. Fractal dimension and FESEM imaging provided evidence for the presence of specific ion effects in BSC aggregate structure and morphology.
The fact that Ca 2+ is more polarized than Mg 2+ under a strong electric field of BSC, implies that Ca 2+ rather than Mg 2+ and K + has a higher probability of appearing in proximity to the charged surface. Therefore, the shielding effect of the electric field is greater in the case of Ca 2+ solution, resulting in a decrease in the activation energy of BSC particle aggregation. Thus, the ion-surface interaction in the electrolyte field should account for cations in the same valence that have significantly different effects on BSC dispersion stability, BSC aggregate structure and surface topography.

CCCs of monovalent and divalent cations in BSC aggregation. According to the DLS experiments,
we made the data into a dynamic view during aggregation and disaggregation, which revealed the interactions and conduced to assessment of velocities of mineral assemblages. The data presented in Fig. 5 show that all of the TAA rates of BSCs in different electrolyte solutions increased significantly at the low concentration ranges but either plateaued or increased gradually at higher concentrations. In each case, the TAA rates of each electrolyte solution can be mathematically expressed by two linear functions both for the low and high concentrations, the two straight lines of the electrolyte concentration will intersect at the intersection point, and it corresponds to CCC. Where the TAA rate was nearly consistent or slightly changed above the CCC, in which the colloid aggregation was DLCA. The portion at which the TAA rate increased linearly with increase in time was RLCA, in which only partial colloid collision resulted in aggregation due to the energy barriers originating from the repulsive force.
The CCC values are found to be equal to 134.3, 13.3 and 4.2 mmol L −1 for KCl, MgCl 2 and CaCl 2 , respectively. The CCC can quantify the ability of different cations to cause soil colloid aggregation. Through analysis of the data of CCCs at various electrolyte concentrations drew a conclusion that CCC changed along with the ionic valence and ion specificity. For coagulation, divalent ions was obviously superior to that of monovalent ions. However, for the same valence cations, CCC Mg was 3.2 times as much as that for CCC Ca . Similar results have been reported in a previous study on mineral montmorillonite aggregation 34 . The CCCs of Ca 2+ and Mg 2+ were 2.38 and 7.99 mmol L −1 , respectively, on montmorillonite aggregation, and CCC Mg was nearly 3 times as much as that for CCC Ca , indicating that the CCC of these two cations in the natural soil particle aggregation process is larger than that in pure montmorillonite because of the presence of organic components in natural soil. The organic components cause a stronger electric field around the soil colloid particles. However, the difference between the natural soil and the mineral particles for CCC Mg is always three times that for CCC Ca , implying that the coagulation ability of Ca 2+ and Mg 2+ is weakened to the same extent. Therefore, the difference in coagulation ability between Ca 2+ and Mg 2+ is not based on the composition of the particle system. Further, it is not clear that whether the difference was caused by the variations in particle size including ionic radius and ion hydration radius? Yet, this difference is contrary to the variation in ion volume reported in previous studies, and the fact that the hydrated radii of Mg 2+ and Ca 2+ were 0.30-0.47 nm and 0.41-0.42 nm, respectively [45][46][47] . Therefore, large differences in CCC for aggregation among this pair of cations did not result from the minute differences in volume because their hydrated radii are also similar. As what the TAA rates and CCC values had indicated, specific ion effects are apparent for different base cations (K + , Mg 2+ and Ca 2+ ). Therefore. For the purpose of indepth understanding of the specific ion effects in the BSC aggregation process, it is thus evaluating the activation energies in KCl, MgCl 2 and CaCl 2 solutions that counts. The functions between the activation energy ∆E(f 0 ) and cation concentration f 0 given above are also plotted in Fig. 6. These data show that no matter what the value of the cation concentration (≤ CCC) is given, the activation   www.nature.com/scientificreports/ specific ion effects change in the order of Ca 2+ < Mg 2+ < K + , and CCC values has a good consistency with the above results on the aggregation rates. By analyzing the activation energies, we can understand intuitively the change trends of the aggregation rates and CCC values in different cation solutions. The K + results show the lowest aggregation rate and the highest CCC value for the aggregation of BSC particles, and the reason for the results is that K + produces the highest activation energies for this aggregation process. As the activation energies decrease from K + to Mg 2+ and finally to Ca 2+ , the corresponding aggregation rates show an increasing trend, whereas the CCC values show a decreasing trend.
Further, the difference in the activation energy of particles between Mg 2+ and Ca 2+ in a certain electrolyte concentration is depicted in Fig. 7. It shows that for Mg 2+ and Ca 2+ , their activation energy differences increase sharply with the decrease of ion concentrations. That is the decrease of ion concentration leads to an obvious increase of specific ion effects. Furthermore, a decrease in electrolyte concentration enhances the electric field near negatively-charged soil particle surfaces. Such a strong electric field enhances cation polarization, and is indicative of the significant role played by the polarization effect. The polarization, also observed in montmorillonite interactions and ion adsorption in previous studies 16,27,48 , indicates that for a real soil system, there is still a nonclassical polarization of cations directly related to the electric field near the soil colloidal surface. Therefore, the ionic polarization at the interface is a non-negligible non-DLVO factor.
The influence of electrolyte concentration on the zeta potential of BSC aggregates. Change in the zeta potential may indirectly reflect the change in the surface charge density. The zeta potential of BSCs were all negative for KCl, MgCl 2 and CaCl 2 solutions within the electrolyte concentration ranges; furthermore, divalent cations in comparison to that of monovalent cations make the zeta potential of BSCs less negative (Fig. 8). Unlike monovalent ions of K + influencing zeta potential primarily through the electric screening effect, divalent ions of Mg 2+ and Ca 2+ adsorb to the surface sites not only by the attractive electrostatic forces but also via the formation of surface complexes with functional groups 49 . However, the zeta potential of BSCs in Mg 2+ was higher than that in Ca 2+ for most electrolyte concentrations, although they have the same ionic valence. Further, the zeta potential is the shear plane potential of the charged colloid, and the distance of the shear plane from the surface in a monovalent cation solution is much greater than that in a divalent cation solution at the same zeta potential 50 . Several studies show that zeta potential is approximately 3.1-6.0 times lower than the surface potential (potential at the original plane of diffusion layer) 50 . The effect of ions on particle surface properties is more easily reflected on the surface potential rather than on the zeta potential.
The data in Fig. 8 shows that at relatively low electrolyte concentration (50 mmol L −1 ) of K + , zeta potential of BSCs is ~ − 27 mV. At this point, the BSC system is in a dispersed stable state, and the stability of the suspension is caused by the repulsive electrostatic forces between colloidal particles. With the increase of cation concentrations, the zeta potential of BSCs became less negative and causes a decreased repulsive force between the particles, indicating unstable suspensions, in agreement with the DLS test results. Colloidal particles, in suspensions with zeta potentials of more than + 20 mV or less than − 20 mV, are stable systems 51 . Some scholars define this threshold value at 25 mV and called it 25 mV rule 52 . Later, researchers estimat the CCCs using 25 mV as a threshold of stability, and they found that 25 mV rule is clearly unreliable for low CCCs (< 100 mV), especially for suspensions containing multivalent ions 18,53 . Here our current results further confirm its limitation.
DLVO simulation of BSC aggregation in K + , Mg 2+ and Ca 2+ electrolytes. The classic DLVO theory, can be used to reveal the stability of colloids, which is affected by van der Waals force and electrostatic force Waals and repulsive electrostatic interactions. In the present study, the CCC values of K + , Mg 2+ and Ca 2+ were obtained, and calculate the sum of DLVO pressure when f 0 = CCC, and plot the result in Fig. 9. When f 0 = CCC, www.nature.com/scientificreports/ the effective collision probability of the particles equals one; that is, aggregation will occur once the particles collide. Therefore, the repulsive barrier between the particles is zero at that time. However, as shown in Fig. 9, the energy barrier ΔW is obviously not zero at CCC K , CCC Ca and CCC Mg according to the classical DLVO theory. Meanwhile, there will be a certain deviation between the aggregation results and the classic DLVO simulation. This deviation value can verify that non-DLVO factors also play a key role in the aggregation process. While it may be difficult to accurately analyze what the non-DLVO factors are, the impact of non-DLVO factors has been reported by others 20 . While considering the non-DLVO factors, it is also necessary to consider the additional energy in the ionsurface interaction (Eq. 11) where w i (total) is the total ion-surface interaction energy of the ith ion species, Z i Fφ 0 is the classic electrostatic interaction energy, w i (add) represents the additional interaction energy of nonclassical DLVO factors, and γ i is the effective ionic charge coefficient of the ith ion species. Regardless of which forementioned effects cause the additional interaction energy of the cation to be changed, the contribution of the total ion-surface interaction energy w i (total) can be measured by the change in the number of charges. Thus, the cationic apparent charge changes from Z i into γ i Z i in the modified calculation. In the aggregation process of a given charged particle, the van der Waals attraction P vdW is a certain value, and the reliable estimation of the energy barrier ΔW relays on the correct electrostatic repulsion, P EDL , regardless of the cation. Within Eqs. (4)- (8), Z i is replaced by γ i Z i to describe the modified electrostatic repulsive interaction. The γi is a key parameter that can determine the  www.nature.com/scientificreports/ strength of additional interaction energy of the ith ion species, and further determines the energy barrier for particle aggregation.
With a series of assumed γ i values, the modified net DLVO force pressure, P DLVO (λ), in CCC K , CCC Ca and CCC Mg are quantitatively calculated with a combination of Eqs. (3)- (10). Wherein, the coefficient makes the modified energy barrier equal to zero, ΔW = 0, is sought after from the γ K , γ Ca , and γ Mg by taking the non-DLVO factors into account. The data in Fig. 10 show that the energy barrier of zero is derived using γ i of 1.31, 1.10, and 1.39 for BSC aggregation in K + , Mg 2+ and Ca 2+ electrolytes, respectively, as non-DLVO factors are considered.
The charge numbers considered non-DLVO factors were in the order γ Ca Z Ca (2.78) > γ Mg Z Mg (2.20) > γ K Z K (1.31) in natural BSC suspension. In addition, the same cation has different charge coefficients in different systems. In the pure mineral model system, the charge numbers of Ca 2+ are 3.38 and 3.92 in montmorillonite and humic acid-montmorillonite complexes and that of Mg 2+ is 2.38 in montmorillonite, respectively 16,34 . The difference of charge numbers between the model system and natural soil system may due to two reasons. Firstly, the colloid surface types are different. The surface of montmorillonite is dominated by permanent charges. While, the surface of natural black soil is more complex including both permanently and variably charged surface. Secondly, the composition of colloidal substances is different. BSC contains clay minerals such as montmorillonite and illite, as well as organic matter, metal oxides and other substances. Studies have shown that divalent ions have different affinities with the surface of humic acid, leading to different aggregation process 19 . Researchers proposed that there is a similar non-classical polarization of anions in the yellow soil colloid system 27 . They concluded that interaction of anions with cations and proton affinities in the electric field was the influencing factors for anion specificities. Therefore, it is necessary to carry out research on soil colloid aggregation with various properties. This means that the effects of non-DLVO factors vary in different systems. All the non-DLVO factors combine to dominate the colloid particle behavior. The effective ionic charge coefficient fitting based on DLS determination accurately quantifies the intensity of non-DLVO forces during natural soil particle aggregation. This calculation provides a basis for further quantitative analysis of non-DLVO factor effects in other processes. The effective charge number described above can be used to supplement the calculation of related parameters considering the non-DLVO effects in the same system.

Conclusions
K + , Mg 2+ and Ca 2+ cations strongly affect the BSC dispersion stability. Ion-surface interactions in a strong electrolyte field is a primary mechanism of cation effects on soil colloid interactions, causing the difference in colloid interaction energy and further affecting soil colloid dispersion stability. The data presented support the conclusion that: (1) pronounced specific ion effects were observed in BSC aggregation follow the sequence of Ca 2+ > Mg 2+ > K + . This sequence is in accord with the activation energy and zeta potential for BSC aggregation in different electrolytes; (2) the effective ionic charge numbers 1.31, 2.20, and 2.78 obtained from the ion specific particles' aggregation kinetics in the presence of K + , Mg 2+ and Ca 2+ can be used to quantify all the non-DLVO factors effect on natural soil colloid aggregation; (3) the specific ion effects of cations further lead to the difference in the structure of the aggregates by affecting the aggregation process, and final affecting the morphology of the aggregates. The new background mechanism above strengthen the understanding of the processes affecting improvement of the stability of black soil aggregates.